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ABSTRACT 

By analysing the Chandra and XMM-Newton archived data of the nearby galaxy cluster Abell 
3158, which was reported to possess a relatively regular, relaxed morphology in the X-ray 
band in previous works, we identify a bow edge-shaped discontinuity in the X-ray surface 
brightness distribution at about 120hj^ kpc west of the X-ray peak. This feature is found to 
be associated with a massive, off-centre cool gas clump, and actually forms the west boundary 
of the cool clump. By calculating the thermal gas pressures in the cool clump and in the free- 
stream region, we determine that the cool gas clump is moving at a subsonic velocity of 
700^,340 ^ O.Gj^Q g) toward west on the sky plane. We exclude the possibility 

that this cool clump was formed by local inhomogeneous radiative cooling in the intra-cluster 
medium, due to the effectiveness of the thermal conduction on the time-scale of ^ 0.3 Gyr. 
Since no evidence for central AGN activity has been found in Abell 3158, and this cool clump 
bears many similarities to the off-centre cool gas clumps detected in other merging clusters in 
terms of their mass, size, location, and thermal properties (e.g. lower temperature and higher 
abundance as compared with the environment), we speculate that the cool clump in Abell 
3158 was caused by a merger event, and is the remnant of the original central cool-core of 
the main cluster or the infalling sub-cluster. This idea is supported not only by the study of 
line-of-sight velocity distribution of the cluster member galaxies, but also by the study of gas 
entropy-temperature correlation. This example shows that the appearance of such massive, 
off-centre cool gas clumps can be used to diagnose the dynamical state of a cluster, especially 
when prominent shocks and cold fronts are absent. 

Key words: galaxies: clusters: individual (Abell 3158) - X-rays: galaxies: clusters- inter- 
galactic medium - kinematics and dynamics 



1 INTRODUCTION 

In the frame of hierarchical clustering scenario, galaxy clusters 
grow in size by merging with subunits, with each major merger 
event lasting for 2 — 5 Gyr (e.g. Roettiger, Loken & Burns 1997; 
Ascasibar & Markevitch 2006; Poole et al. 2006). In such major 
mergers both ram pressure stripping and slingshot can generate re- 
markable X-ray substructures in gas density and temperature distri- 
butions, which are often accompanied by shocks and/or cold fronts 
that exhibit arc-shaped or edge-like morphologies (e.g. Markevitch 
& Vikhlinin 2007 and references therein). These substructures are 
expected to contain valuable information not only on merger dy- 
namics itself, but also on thermal and chemical evolutions of the 



hot intra-cluster medium (ICM), which helps understand the nature 
of dark matter (Markevitch et al. 2004; Randall et al. 2008). 

Theoretical studies and numerical N-body simulations (e.g. 
Roettiger et al. 1997; Ricker 1998; Poole et al. 2006) showed that 
in a typical major merger event temperature substructures may sur- 
vive for a relatively long time of fthcrmai — ngi^skT/Ug^A, where 
rigas, T, and A are gas density, temperature, and cooling func- 
tion, respectively. The derived tthermai (up to about 4 Gyr in some 
cases) is much longer than the lifetimes of shocks and cold fronts, 
since shocks and cold fronts are usually smeared out quickly on 
dynamic time-scales characterized by the sound crossing time. Ac- 
tually, high-quality Chandra and XMM-Newton observations have 
shown that in merging clusters there is a prevailing existence of 
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massive, off-centre cool gas clumps, which are about 50 — 300 kpc 
in size, reside within the central (<600 kpc) region of the cluster, 
and possess an average temperature 1 — 4 keV lower than the am- 
bient gas. 

In IE 0657 - 56 (the Bullet cluster; Markevitch et al. 2002), 
Abell 520 (Govoni et al. 2004), Abell 2256 (Sun et al. 2002), and 
Abell 3667 (Vikhlinin, Markevitch & Murray 2001), the off-centre 
cool gas clump is found to be moving adiabatically behind a bow 
shock and/or a cold front, and thus can be directly identified as 
the remnant of the original central cool-core of the infalling sub- 
cluster. On the other hand, in clusters such as Abell 754 (Marke- 
vitch et al. 2003), Abell 2065 (Chatzikos, Sarazin & Kempner 
2006), and Abell 2255 (Sakelliou & Ponman 2006), the off-centre 
cool gas clump is not escorted by either a shock or a cold front. 
However, unambiguous evidence obtained in radio, optical, and X- 
ray bands still indicates that the cool gas clump was driven out of 
the core region of either the main cluster (Abell 754; Markevitch 
et al. 2003; Henry, Finoguenov & Briel 2004) or the infalling sub- 
cluster by the merger. 

In fact, the off-centre cool gas clumps in above two types of 
clusters share many similarities in terms of their sizes, locations, 
and gas masses (~ 10^^ Mq). Therefore, an interesting ques- 
tion arises: can the appearance of such massive, off-centre cool 
gas clumps be safely used to diagnose the dynamical state of a 
cluster, even when there is a lack of violent merger signatures, 
such as prominent shocks, cold fronts, and other associated ra- 
dio/optical substructures? Clearly, in order to address this issue, we 
need to examine gas temperature distribution in more galaxy clus- 
ters that show only weak merger signatures and have evolved into 
late merger stages (usually mean post core passages). 

In this work we present a detailed joint Chandra and XMM- 
Newton study of the nearby rich galaxy cluster Abell 3158 (rich- 
ness class 2, Quintana & Halven 1979; z — 0.0597, Struble & 
Rood 1999), which contains 105 member galaxies within the cen- 
tral 2.4/17/ '^P'^ (Biviano et al. 2002). In literature (e.g. Irwin, 
Bregman & Evrard 1999; Lokas et al. 2006; Chen et al. 2007), this 
cluster is usually referred as a relaxed, non-cool-core system, be- 
cause neither significant X-ray substructures nor diffuse radio emis- 
sion has been reported in the previous works (Ku et al. 1983; Mohr, 
Mathiesen & Evrard 1999; Schuecker et al. 2001; Mauch et al. 
2003). In the optical band, however, the cluster shows some unusual 
properties that can be possibly related to a recent merger event, as 
are summarized as follows. First, the cluster hosts three luminous 
elliptical galaxies that are brighter than 9 x 10^" Lb,q (Paturel et 
al. 2003), all classified as a cD galaxy by Sersic (1974). Two of 
them, the brightest member PGC 13641 (EO, Lb = 1.45 x lO" 
Lb,©) and the SO galaxy PGC 13652 (Lb = 9.39 x 10^° Lb,©), 
constitute a closely aligned central galaxy pair, whose projected 
separation is 81.3/if/ kpc, and jointly dominate the cluster (Quin- 
tana & Havlen 1979). The second luminous member PGC 13679 
(SO, Lb = 9.48 x 10^" Lb,©) is located at about 417.8/if/ kpc 
southeast of the central galaxy pair, and dominates a local sub- 
cluster whose recessional velocity is about 1 300 km s^^ larger 
than the mean recessional velocity of the main cluster (17 500 km 
s~^; Lucey et al. 1983; Kolokotronis et al. 2001). Second, as shown 
in figure 7 of Smith et al. (2004), the line-of-sight velocity distri- 
bution of the cluster's member galaxies exhibits a plateau at the 
high-velocity end, which is often seen in merging systems. Third, 
Kolokotronis et al. (2001) reported that in this cluster there is an 
offset of about 200/if/ kpc between the optical flux-weighted cen- 
troid and the optical density peak. Based on these we speculate that 



Abell 3158 is evolving at a late merger stage and is thus an ideal 
target for our investigation. 

Throughout the paper, we adopt the cosmological parameters 
//o = 71 km s"^ Mpc"\ Qm = 0.27, and = 0.73, so that 
1" corresponds to about 1.14/i^/ kpc at the redshift of the cluster. 
We utilize the solar abundance standards given by Grevesse and 
Sauval (1998), where the iron abundance relative to hydrogen is 
3.16 x 10^'' in number. Unless stated otherwise, the quoted errors 
are the 90% confidence limits. 



2 OBSERVATIONS AND DATA REDUCTIONS 

2.1 Chandra 

Except for an extremely short exposure on 2007 September 16 (5.1 
ks, ObsID 7688), which is not used in this work, Abell 3158 has 
also been observed with Chandra on 2002 June 19 (31.4 ks, Ob- 
sID 3712) and June 21 (25.1 ks, ObsID 3201), respectively, with 
chips 0, 1,2, 3, and 6 of the Advanced CCD Imaging Spectrom- 
eter (ACIS) operating in VFAINT mode. We used the Chandra 
data analysis package CIAO version 4. 1 and C ALDB version 4. 1 .2 to 
process the archived data in the standard way, by starting with the 
level- 1 event files. Corrections for the charge transfer inefficiency 
and time dependent gain have been applied. We kept events with 
ASCA grades 0, 2, 3, 4, and 6, and removed all the bad pixels, bad 
columns, and columns adjacent to bad columns and node bound- 
aries. We examined the 0.3 — 10.0 keV lightcurves extracted from 
the background regions defined on ACIS chips and 1 , and found 
that there are no strong flares that increase the background count 
rate to > 120% of the mean quiescent value. The obtained net ex- 
posures are 30.9 ks and 24.8 ks for the two observations, respec- 
tively. In the spectral analysis, we extracted the Chandra spectra in 
0.7 — 8.0 keV from the two observations separately, and fitted them 
simultaneously with the same model, except that their normaliza- 
tions were left free. Background spectra were extracted from the 
Chandra blank-sky fields; a crosscheck based on the use of local 
background yielded essentially the same results. 

2.2 XMM-Newton 

Abell 3158 was observed with XMM-Newton on 2005 November 
22 (22.4 ks, ObsID 0300210201) and 2006 January 18 (9.4 ks, Ob- 
sID 030021 1301), respectively, with the European Photon Imaging 
Camera (EPIC) operating in PrimeFullWindow (MOS 1 and M0S2) 
and PrimeFuUWindowExtended (pn) modes. Although during the 
second observation the Reflection Grating Spectrometer (RGS) was 
also turned on, in this work we limited our analysis to the first ob- 
servation only by using the latest version of SAS (version 8.0.1) and 
its standard filterings. We kept the events with FLAG=0 and PAT- 
TERNS 0-12 for MOSl and M0S2, and the events with FLAG=0 
and PATTERNS — 4 for pn. We removed all the time intervals con- 
taminated by soft proton flares, during which the 10 — 12 keV count 
rate exceeds the 2(t limit of the mean quiescent value (see, e.g. 
Katayama et al. 2004). The cleaned MOSl, M0S2, and pn datasets 
have effective exposure times of 20.9 ks, 20.6 ks, and 12.2 ks, re- 
spectively. In the spectral analysis, we removed all the bright X-ray 
point sources and adopted conservative energy cuts at 0.5 keV and 
8.0 keV. The obtained MOSl, M0S2, and pn spectra were fitted 
simultaneously using the same model, except that their normaliza- 
tions were left free. Backgrounds spectra were generated from the 
blank-sky event lists, which were filtered in advance using the same 
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selection criteria (PATTERN, FLAG, etc.) as used for the source spec- kpc and thus the edges have been confirmed as cold fronts caused 
tra. by mergers (Markevitch et al. 2002; Vikhlinin et al. 2001). 



3 X-RAY ANALYSIS AND RESULTS 

3.1 X-ray surface brightness discontinuity 

In Fig. la and lb, we show the Chandra ACIS and combined 
XMM-Newton EPIC images of the central Mpcxl.5/ifi^ 
Mpc (21'. 9 X 21'. 9) region of Abell 3158 in 0.3 - 2.0 keV, re- 
spectively, which have been adaptively smoothed and corrected 
for both exposure and vignetting. The XMM-Newton image was 
generated by combining the MOSl, M0S2, and PN events to- 
gether using the SAS 'image' script. The X-ray peak of the clus- 
ter (R.A.=03h42m52.7s Dec=-53d37m37.2s J2000) is found at 
about n.Shjl kpc (0'.26) northwest of the optical centroid of 
PGC 13641 (R.A.= 03h42m52.9s Dec=-53d37m52.2s; Katgert et 
al. 1998), the brighter member of the central galaxy pair. The corre- 
sponding optical Digital Sky Survey (DSS) image is shown in Fig. 
Ic, on which the Chandra intensity contours obtained from Fig. la 
are overlaid. 

On both Chandra and XMM-Newton images it can be clearly 
seen that the X-ray isophotes are elongated in nearly the east-west 
direction with an ellipticity of ~ 0.3. And on the Chandra image a 
bow edge can be identified at about 120h^^ kpc (l'.88) west of the 
X-ray peak. According to its shape and location on the detector, the 
possibility of the feature being associated with the artifacts caused 
by CCD gaps or bad columns can be excluded. 

In order to examine the significance of the bow edge in a 
quantitative way, we define two sets of partial elliptical annuli, 
which are confined in the east and west sector regions shown in 
Fig. 2a, and extract the exposure-corrected X-ray surface bright- 
ness profiles (SBPs) with both Chandra ACIS and XMM-Newton 
EPIC (Fig. 26). Because the gas temperature of Abell 3158 is suf- 
ficiently high (5.8 keV; Reiprich & Bohringer 2002), by following, 
e.g. Markevitch et al. (2000), we restricted the SBP extractions in 
0.3 — 2.0 keV to minimize the dependence of X-ray emissivity on 
possible temperature fluctuations (see §3.2). We attempt to apply 
the standard /? model to describe the extracted SBPs, which is ex- 
pressed as 

S{r) = So[l + {r/r,)T'~'^ + ^bkg (D 

(Jones & Forman 1984), where is the core radius, /3 is the slope, 
5'bkg is the background, and Sq is the normalization. We find that 
in the east sector the SBPs can be fitted very well with the /3 model 
(Table 1), whereas in the west sector the SBPs show a significant 
central excess beyond the best-fit /3 model for the outer (> 120h^^ 
kpc) regions. Within about 90/ifi^ - kpc (l'.32 - l'.68), 

where the bow edge is identified visually, there exists a sudden dis- 
continuity, across which the Chandra and XMM-Newton surface 
brightness increases inwards by a factor of ~ 1.6 and ~ 1.3, re- 
spectively; this difference may be caused by the different distri- 
butions of the point spread function (PSF) of Chandra ACIS and 
XMM-Newton EPIC. These results confirm the existence of the bow 
edge shown in Fig. la, and indicate that it is not related to the cen- 
tral surface brightness excess usually seen in the cD galaxies (Mak- 
ishima et al. 2001), which are not associated with a surface bright- 
ness discontinuity. Also, we note that the X-ray surface brightness 
discontinuity in Abell 3158 is similar to, although not as sharp as, 
those found in IE 0657 — 56 and Abell 3667, in which the jumps 
of the surface brightness at the edges are 2 — 3 within ~ lOh^^ 



3.2 An unusual off-centre cool gas clump 

In order to investigate the thermal properties of the bow edge, we 
calculate the two dimensional gas temperature distribution in the 
central 0.6/ifi^ Mpcx0.6/ifi^ Mpc (8'. 77 x 8'. 77) region of AbeU 
3158, by following the approach of, e.g. O'Sullivan et al. (2005) 
and Gu et al. (2009). In the calculation only the XMM-Newton EPIC 
datasets are used, because the data statistics of the Chandra ACIS 
datasets is not good enough for us to map the temperature distribu- 
tion. To be specific, we first define about 4000 cells in the 0.6h^^ 
Mpcx0.6/if/ Mpc region, whose centres (r^, i = 1,2,3...) are 
randomly distributed with a separation of < 10" between the cen- 
tres of any two adjacent cells. Each cell is apportioned with an 
adaptive radius of 0' .3 — 0' .8, so that it encloses > 3000 photons in 
0.5 — 8.0 keV after the point sources are excluded. Then we extract 
the MOS 1 , MOS2, and pn spectra from each cell and fit them simul- 
taneously with an absorbed APEC model coded in XSPEC vl2.4.0. 
The redshift and absorption are fixed to z = 0.0597 and the Galac- 
tic value Nu = 1.62 x 10^° cm"^ (Dickey & Lockman 1990), 
respectively. After the obtained best-fit gas temperature is assigned 
to the corresponding cell (Tc(ri)), and this process is repeated for 
all cells, we calculate the projected temperature at any position r as 

r(r) =^[Gr,(i?r.rJrc(r,)]/X]Gr,(i?r.rJ, Rr,r, <s(r,) 

r J 

(2) 

(Gu et al. 2009), where Rr,Ti is the distance between r and the cen- 
tre of the cell i, s{ri) is defined as the radius of cell i, and Gn is the 
Gaussian kernel with a scale parameter a = s(ri). The obtained 
temperature map, along with the la error map that is calculated in 
nearly the same way as the temperature map, are shown in Fig. 3a 
and 36, respectively. 

We find that, despite the relatively regular appearance of the 
X-ray images (Fig. 1), the spectral distribution of the projected 
gas temperature is highly asymmetric, inferring that the cluster has 
not recovered from a violent event on cluster scales. In particu- 
lar, we find that there exists a cool gas clump at about SAh^-^ kpc 
(l'.23) west of the X-ray peak, whose linear scale is about 80h^^ 
kpc (1'.17). The projected gas temperature within this cool clump 
ranges from about 4.4 keV to 5.0 keV, which is apparently lower 
than that of the ambient gas (5.2 — 6.0 keV). Most interestingly, 
the west boundary of the cool gas clump coincides with the bow 
edge perfectly, indicating that they are likely to have the same ori- 
gin. 

To clarify the significance of the existence of the cool gas 
clump, we define four sector regions in Fig. 3a, two of which have 
the same angles as those defined earlier in Fig. 2a, and divide each 
of these sector regions into a series of wider partial elliptical an- 
nuli. We extract both the Chandra ACIS and XMM-Newton MOSl, 
M0S2, and pn spectra from these annuli, and fit them with an 
absorbed APEC model when the redshift and absorption are fixed 
again as above; allowing the absorption to vary does not improve 
the fits. In the fittings, the Chandra ACIS spectra extracted from the 
two observations are treated as group one, and the XMM-Newton 
MOSl, M0S2 and pn spectra are treated as group two. Each group 
of spectra are fitted simultaneously with the same model param- 
eters, except for that the normalizations are left free. The derived 
projected temperature profiles in the four sectors are plotted in Fig. 
4a — 4d. It can be clearly seen that the Chandra and XMM-Newton 
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temperatures are consistent with eacli other (68% confidence level), 
and in the west sector the temperature in 45.4 — 95.3h^^ kpc 
(4.85 ± 0.17 keV), where the cool gas clump is located, is lower 
than those of the adjacent regions at 68% confidence level (see also 
Table 2); no such temperature drop is found in other three direc- 
tions. This is further confirmed by the results of the deprojected 
analysis of the XMM-Newton spectra for the west sector (Table 2 
and Fig. 4e — 4/), which shows that the gas temperature of the 
cool gas clump (S.OSIq *'^ keV for W2 region) is actually lower 
than those of the innermost region (5.71^0 71 keV for Wl region) 
and the next outer region (5.38!to'5i keV for W3 region) at the 
90% confidence level, and cannot be ascribed to the uncertainties 
in determining the metal abundance. 



3.3 Faint Cold front associated witli the cool gas clump 

To study the nature of the cool gas clump, we carry out a simple 
hydrodynamic study to determine the velocity of the edge in the 
environment. By applying the best-fit deprojected spectral parame- 
ters (Table 2), here we revisit the X-ray surface brightness profiles 
of the west sector (Fig. 2b) and fit it again with two models. In 
the first model (model A), the spatial distribution of gas density is 
described with two j3 components as 



(7?) = {n^,i[l+(i?/i?ci) 



'"''+ng,2[l+(i?/i?ci 



3/321 1/2 



(3) 

where R is the 3-dimensional radius, Rc is the core radius, and /? 
is the slope. In the second model (model B), the gas density distri- 
bution is broken at the truncation radius i?cut, described with one 
truncated power-law component and one truncated jS component as 



ngaa(-R) 



ng,i(i?/i?cut)-" 
n^,2[l + {R/R.f]-^^^^ 



R < Rcut 
R ^ Rcut 



(4) 



Acceptable fits to both the Chandra and XMM-Newton SBPs can 
be obtained with model B only (Table 1). The Chandra best-fit sug- 
gests a density jump by a factor of /jump = 1.7 ± 0.1 at _Rcut ~ 
111.0 ± AAh^^ kpc (roughly the position of the edge), and the 
XMM-Newton best-fit gives consistent values of /jump = 1.5 ±0.1 
and i?cut = 104.9 ± 2.2hrl kpc. 

Following Vikhlinin et al. (2001), we use the best-fit gas den- 
sities across the edge (Eq.(4)) and gas temperatures for W2 and W3 
regions (§3.2 and Fig. 4e) to calculate the thermal gas pressures in 
the cool gas clump (Pq) and in the free-stream region (Pi), respec- 
tively. The former pressure is assumed to be in pressure equilibrium 
with the gas pressure at the stagnation point (denoted as the place 
where the relative gas velocity vanishes). The pressure ratio is es- 
timated to be Po/Pi = 1.3 ± 0.2 (68% confidence level). This 
allows us to determine the Mach number (A^) of the cool clump 
to be O.GIq 3, which corresponds to a velocity of 7OOIJ4Q km s~^. 
This subsonic velocity of the clump, along with the results obtained 
in §3.1 (i.e., surface brightness discontinuity) and §3.2 (i.e., gas 
temperature and density discontinuities), allow us to conclude that 
there exists a faint cold front at the west boundary of the cool gas 
clump. 



4 VELOCITY PLATEAU 

In early works of Biviano et al. (1997 and 2002), Girardi et al. 
(1996), and Kolokotronis et al. (2001), both the two dimensional 
spatial distribution and the line-of-sight velocity distribution of the 
member galaxies in Abell 3158 were studied, and no significant 



substructure was reported. However, after tighter and more reliable 
selection criterion for member galaxies, which was based on the 
galaxy colour-magnitude relation (e.g. Colless & Dunn 1996; Fer- 
rari et al. 2003; Maurogordato et al. 2008), was applied. Smith et 
al. (2004) re-constructed the line-of-sight velocity distribution with 
higher fidelity and found that there exists a visible plateau beyond 
the Gaussian distribution at the high-velocity side. Since Smith et 
al. (2004) did not use a quantitative approach to assess the signifi- 
cance of the high-velocity plateau, in which the merger information 
might be contained, here we draw the velocity data of the member 
galaxies in Abell 3158 from Smith et al. (2004), and re-analyse the 
line-of-sight velocity distribution profile. After the model fitting of 
this distribution is finished, we also attempt to investigate the 2- 
dimensional spatial distribution of the galaxies in the high-velocity 
plateau. 

First, we plot the line-of-sight velocity distribution of the 68 
member galaxies identified by Smith et al. (2004) in Fig. 5a, which 
shows that the high- velocity plateau is located at about 19 000 km 
s~^. We fit the observed distribution with a single Gaussian pro- 
file, and calculate the Kolmogorov— Smimov statistic for the ob- 
served distribution against the best-fit Gaussian model {y^ /dof — 
31.6/9). We find that the observed distribution has a probability 
of < 10% of being Gaussian. We then attempt to fit the observed 
distribution with a two-Gaussian model. The best-fit {X^ /dof — 
9.8/6) gives an average velocity of < vi >— 17 380 ± 50 km s^^ 
and a corresponding variance of cr„,i = 520 ± 50 km s~^ for the 
main Gaussian component, and < V2 >= 19 120 ± 130 km s~^ 
and (j„,2 = 510± 100 km s~^ for the high-velocity plateau. By ap- 
plying the F-test and the Kaye's Mixture Model (KMM; McLach- 
lan & Basford 1988; Ashman, Bird & Zepf 1994) test, the latter of 
which is based on a maximum likelihood algorithm, we find that 
the second Gaussian component is required at 94.3% confidence 
level and preferred at a significant probability of 90%, respectively. 

To investigate if the galaxies in the high-velocity plateau form 
a real substructure in the cluster, we divide the member galaxies 
into two subgroups: one low-velocity (15 300 — 18 500 km s~^) 
subgroup and one high- velocity (18 500 - 20 000 km s"^) sub- 
group, which consist of 51 and 17 galaxies, respectively. According 
to the best-fit two-Gaussian model, these two subgroups roughly 
corresponds to the two Gaussian components, respectively, with up 
to about two of the galaxies in the high- velocity subgroup coming 
from the main Gaussian component. We find that the galaxies be- 
longing to the high- velocity subgroup are distributed mostly in the 
central SOOh^l kpc (7'. 3), as well as the west part of the cluster, 
with a geometric centroid of the high-velocity galaxies is located 
at about SSbh^i kpc (5'. 2) northwest of the X-ray peak, which is 
not coincident with any X-ray substructures or luminous member 
galaxies. The high-velocity subgroup roughly includes the member 
galaxises of the PGC 13679 sub-cluster, since the member galaxies 
and the boundary of the latter are difficult to be determined by 2- 
dimensional spatial distribution alone (e.g. Colless & Dunn 1996; 
Boschin et al. 2006). The galaxies in the low-velocity subgroup that 
is dominated by the central PGC 13641 -PGC 13652 galaxy pair, on 
the other hand, are scattered symmetrically in the field. These re- 
sults suggest that the galaxy velocity separation (Fig. 5a; see also 
Smith et al. 2004) has a dynamical nature, which is most likely to 
be a major merger, with a merger mass ratio of about 1 : 1 — 3, as 
estimated from the velocity variance ratio (cr„^2/o"i>,i — 1:1) and 
the galaxy number ratio (~ 17/51 = 1 : 3) of the two subgroups. 

In literature Abell 3158 was classified as a 'single-component' 
and 'relaxed' system (Kolokotronis et al. 2001; Lokas et al. 2006), 
which possesses a relatively regular morphology in the X-ray band 
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(Ku et al. 1983; Mohr et al. 1999; Schuecker et al. 2001). How- 
ever, considering none of violent merger signatures is detected in 
the cluster, the existence of the high-velocity subgroup indicates 
that the cluster is evolving into the late stage of a major merger 
event. This may explain the origin of the off-centre cool gas clump 
detected in the west sector region, which is moving at a subsonic 
velocity (§3.2— §3.3; see also §5.1). 



5 DISCUSSION 

5.1 Merging origion of the massive, off-centre cool gas clump 

The most prominent X-ray feature detected in Abell 3158 is the 
massive, off-centre cool gas clump, which is moving behind a faint 
cold front. It does not resemble the cool substructures caused by 
AGN activity (e.g. Forman et al. 2005; Nulsen et al. 2002), which 
are usually tightly associated with the low-density regions such as 
the AGN-induced radio lobes and X-ray cavities, since no evidence 
for AGN activity and its remnants has been reported in Abell 3158 
according to the existing multi-band observations. In such a sense, 
we speculate that the cool clump was formed either by the inho- 
mogeneous radiative cooling in the ICM, or by the ram pressure 
stripping or slingshot during a merger event. The first possibility 
can be ruled out immediately, because given the best-fit spectral 
parameters for the cool clump, the radiative cooling time of the gas 
in the clump (cf. Sarazin 1986) is 



3.5 X 10"^ yr x 



2 X 10-» cm-3 



kT 
lOkeV 



1/2 



9.8 Gyr, 



(5) 



where continuum and line emissions are both included. This is 
close to the Hubble time (~ 12.9 Gyr) of the cluster, and is sig- 
nificantly longer than the Coulomb conduction time 



tr. 



9.2 X 10" yr x 



100 kpc 



2 X 10-3 cm-3 

5/2 



lOkeV 



0.06 Gyr (6) 



(Markevitch et al. 2003), where / — 80/if/ kpc is the linear scale 
of the cool clump, fcTcentre = 5.5±0.1 keV is the average gas tem- 
perature for the cluster's central 200h^^ kpc (2'. 92) region. If the 
conduction is suppressed by a factor of 5 due to the tangled mag- 
netic field in the ICM (e.g. Narayan & Medvedev 2001; Zakamska 
& Narayan 2003; Chandran & Maron 2004), the effective conduc- 
tion time increases to 0.3 Gyr. Consequently, even if the cool gas 
clump was formed through radiative cooling, it would have been 
destroyed very quickly by conduction, which indicates that the cool 
gas clump is most likely caused by merger. The same conclusions 
also hold for off-centre cool gas clumps found in other merging 
clusters (Table 3). 

Assuming that the cool gas clump in Abell 3158 has a spher- 
ical geometry with a diameter of 80 — IWh^^ kpc (Fig. 3a), 
and has a constant density, the gas mass is estimated to be about 
3.0 - 7.7 X 10^" Mq. This value is comparable to masses of the 
off-centre cool clumps identified in other merging clusters, which 
range from about 2.0 x 10^° M@ to 1.3 x 10^^ Mq (Table 3). It also 
falls into the gas mass range of the cool-core regions of poor clus- 
ters and massive groups (e.g. Mulchaey et al. 1996; Gastaldello et 
al. 2007; Baldi et al. 2009). In fact, as shown in Table 3, not only the 



gas mass, but also the size, temperature gradient, and the relatively 
high metal abundance of the cool clump in Abell 3158 resemble 
those of other off-centre cool clumps detected in IE 0657 — 56 
(Markevitch et al. 2002), Abell 754 (Markevitch et al. 2003), Abell 
2065 (Chatzikos et al. 2006), Abell 2255 (Sakelliou & Ponman 
2006), Abell 2256 (Sun et al. 2002), and Abell 3667 (Vikhlinin 
et al. 2001). Note that most of the off-centre cool clumps in these 
merging systems are considered to be the remnant of the cool-core 
of either the main cluster or the infalling sub-cluster. Based on all 
these similarities, it is natural to speculate that the cool gas clump in 
Abell 3158 was formed by the ram pressure stripping or slingshot 
during a merger event. 

The above speculation is enhanced by the study of the cor- 
relation between gas temperature and entropy, which is defined as 
5* — kT/rLgie in terms of gas temperature kT and density rigas. By 
quoting the data of Zhang et al. (2007) and Sanderson, O' Sullivan 
& Ponman (2009), we re-calculate the redshift-corrected gas en- 
tropy at O.I-R500 (^500 is the radius at which the mean over-density 
is 500 with respect to the critical density of the universe) for a sam- 
ple of 13 cool-core clusters and a sample of 19 non-cool-core clus- 
ters. We also calculate the gas entropy of the seven cool gas clumps 
listed in Table 3. The entropy-temperature distribution is plotted 
in Fig. 6, where the strong correlations for the cool-core and non- 
cool-core clusters are characterized by two power-law models, re- 
spectively (Sanderson et al. 2009). It can be clearly seen that cool 
gas clumps in Abell 3158 and most of other merging clusters fol- 
low the relation for cool-core clusters, which reveals their nature as 
cool-core remnants. Note that the gas entropy of the cool clump in 
IE 0657-56 is significantly below the relation for cool-core clus- 
ters at least la limit (Fig. 6), which is possibly due to that the cool 
clump is compressed extremely as it experiences the ram pressure 
stripping force and passes through the deep potential of the main 
cluster during this drastic high- velocity merger event (Markevitch 
et al. 2002). 

Based on above data analysis, calculations, and discussions, 
we propose that the appearance of the massive off-centre cool gas 
clumps can be regarded as a robust diagnostic of merger state. Com- 
bining the study of cold fronts (e.g. Owers et al. 2009), which sur- 
vive longer than the shocks, with the study of such cool gas clumps, 
we may be able to place tight constraints on the merger history. 



5.2 Possible merging scenarios 

In Abell 3158, the X-ray bow edge of the faint cold front is point- 
ing towards west in projection (§3.1 — §3.3), which indicates that 
the symmetrical axis of the merging process is reasonably close to 
the west-east direction, although the projected effect lowers the de- 
tectability of the edge. The direction is also roughly consistent with 
the extending direction of the three cD galaxies (Fig. la), as well 
as the galaxy distribution of the high-velocity subgroup (Fig. 56). 
Straightforwardly, the high- velocity subgroup has already made its 
closest approach to the core of the main cluster roughly from east 
to west. In this scenario, the apparent discrepancy between the 
moving velocity of the cool gas clump (700^340 km s~^; §3.3) 
and the relative line-of-sight velocity of the high- velocity subgroup 
(1740±180 km s~^; §4) demonstrates the gas component of the 
high-velocity subgroup possibly has been significantly slowed by 
the ram pressure stripping force. Regarding the late stage of the 
major merger (§4), perhaps the high-velocity subgroup is in its sec- 
ond core passage. On the other hand, because the three cD galaxies 
all lie well behind the cold front (Fig. la), perhaps the cluster is 
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in another merging scenario, in wliicii the cool gas clump has been 
slingshot out from the centre of the main cluster, as expected af- 
ter core passage, although the optical centroid of the high-velocity 
subgroup is not behind the faint cold front edge. 



6 SUMMARY 

By analysing the Chandra and XMM-Newton data of the nearby 
galaxy cluster Abell 3158, we identify a massive, off-centre cool 
gas clump, which is moving at a velocity of = O.GIq ^, toward 
west on the sky plane behind a faint cold front. The possibilities of 
this substructure being formed by inhomogeneous radiative cool- 
ing and by central AGN activity are excluded. Based on the results 
obtained in the X-ray and optical analysis, we speculate that the 
cool gas clump is the remnant of the original central cool-core of 
the main cluster or the infalling sub-cluster during a major merger 
event. This case shows that the appearance of such massive, off- 
centre cool gas clumps can be used to diagnose the dynamical state 
of a cluster, especially when violent merger signatures (e.g. promi- 
nent shocks and cold fronts) are absent. 
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B 


Chandra 


4.51 ± 0.09 


0.07 ± 0.04 


111. 0± 4.4 


3.16 ±0.05 


0.63 ±0.01 


238.9 ±3.3 


20.7/15 


B 


XMM 


4.72 ± 0.08 


0.03 ±0.03 


104.9 ± 2.2 


3.92 ±0.05 


0.66 ±0.01 


200.0 ±2.1 


33.5/24 



Table 1. Best-fits to the exposure-corrected surface brightness profiles (Fig. 26), which are extracted in 0.3 — 2.0 keV from the partial elliptical annuli defined 
in the east and west sector regions (Fig. 2a). " A single /3 model is sufficient in fitting the SBPs extracted in the east sector (§3.1). Best-fit deprojected 
spectral parameters obtained in §3.2 (see also Table 2 and Fig. 4e) and two gas density models (§3.3; model A=/3 component-H/9 component (Eq. 3); model 
B=tmncated power-law component-l-truncated /3 component (Eq. 4)) are used to fit the Chandra and XMM-Newton profiles extracted in the west sector. 
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Table 2. Gas temperature and abundance distributions in the west sector region. " We divide the west sector region into wider partial elliptical annuli, 
i.e., Wl— W5 regions (Fig. 3a). Equivalent inner and outer radii for each partial elliptical annulus, which are defined as r = \/ab, where a and b are the 
corresponding semi-major and semi-minor axes, respectively. An absorbed APEC model is used to fit the Chandra ACIS and XMM-Newton MOS 1 +MOS2+pn 
spectra extracted in the Wl— W5 regions. The redshift and absorption are fixed to ^ = 0.0597 and the Galactic value Nn = 1.62 X 10^" cm"^ (Dickey & 



Lockman 1990), respectively. EiTors are quoted at the 68% confidence. 
MOSl+MOS2+pn spectra. Errors are quoted at the 90% confidence. 



The PROJCT model coded in XSPEC vl2.4.0 is used to deproject the XMM-Newton 
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Table 3. Properties of the massive, off-centre cool gas clumps, which are detected in IE 0657-56 (Markevitch et al. 2002), Abell 754 (Markevitch et al. 2003; 
Henry et al. 2004), Abell 2065 (Chatzikos et al. 2006), Abell 2255 (Sakelliou & Ponman 2006), Abell 2256 (Sun et al. 2002), and Abell 3667 (Vikhlinin et al. 
2001; Briel et al. 2004), respectively. Cooling times of the cool clumps are calculated from Eq (5), and effective conduction times are esitmated from Eq (6) 
by assuming a tangled magnetic effect. Gas masses of the cool clumps are calculated by assuming that the gas density, temperature, and metal abundance are 
constant in the clump. Gas entropies are given by 5 = kT/rigis. 
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Figure 1. (a)-(6) Adaptively smoothed and exposure-corrected Chandra ACIS-I (ObsID 3712) and XMM-Newton EPIC (ObsID 0300210201) images of the 
central \.5h~l Mpcx l.S/if/ Mpc (21'.93 X 21'.93) region ofAbell 3158. Both images are extracted in 0.3—2.0 keV and plotted in logarithmic scale. One 
small box and three crosses are used to mark the peak of the X-ray halo and the centres of luminous member galaxies PGC 13641, PGC 13652, and PGC 
13679, respectively, (c) Optical DSS B-band image for the same sky field as (a) and (b), on which the Chandra X-ray intensity contours (square root scale) 
are plotted. 
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Figure 2. (a) East and west sector regions used to extract the X-ray SBPs, which are shown in (6). Here the Chandra ACIS image is used as the background 
image. (6) Chandra ACIS and XMM-Newton EPIC 0.3—2.0 keV SBPs extracted from the partial elliptical annuli defined in the two sector regions. For SBPs 
extracted in the east sector, a single /3 model can give an acceptable fit to the observed SBPs (dashed lines; §3.1 and Table 1). For SBPs extracted in the west 
sector, dashed, dotted, and solid lines are used to show the best-fits obtained with the single /3 model (§3.1), two-/3 model, and truncated powerlaw-H/3 model 
(§3.3 and Table 1), respectively. 
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Figure 3. (a) Projected gas temperature map for the central Q.(yh^^ Mpcx0.6/i^j^ Mpc (8'. 77 X 8'. 77) region, along with the 68% error map (b), which are 
calculated from the XMM-Newton MOSl+MOS2+pn data. Partial elliptical annuli defined in the north, east, south, and west sector regions are shown in (a), 
which are used to extract both the Chandra and XMM-Newton spectra to be analysed in §3.2. 
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Figure 4. (a)—(d) Projected gas temperature profiles obtained witii Chandra ACIS (cross) and XMM-Newton MOSl+MOS2+pn (diamond) data for tlie nortli, 
east, soutli, and west sectors defined in Fig. 3a, along with the 68% errors, (e) Deprojected XMM-Newton temperature profiles for the west sector, along with 
the 68% (solid) and 90% (dashed) en'ors (see also Table 2). (/) Fit-statistic contours of temperature and abundance at the 68% (Ax^ = 2.30; solid) and 90% 
(Ax^ = 4.61; dashed) confidence levels for the inner three partial annuli of the west sector, all obtained in the XMM-Newton deprojected analysis. 
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Figure 5. (a) Line-of-sight velocity distribution of 68 member galaxies identified in Abell 3158. Original velocity data are drawn from Smith et al. (2004). The 
distribution is fitted with a two-Gaussian model (solid) and a single Gaussian model (dashed; §4), respectively, (b) DSS optical image for the central 3.0h^-^ 
Mpcxi.OhZ^ Mpc (43'. 86 X 43'. 86) region of Abell 3158, where the low-velocity (v < 18 500 km s~i) and high-velocity {v > 18 500 km s~i) member 
galaxies are marked with small and large circles, respectively. The X-ray peak is marked with a small box and the geometric centroid of the high-velocity 
subgroup is marked with a cross. 
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Figure 6. Redshift-corrected gas entropy measured at O.liJsoo (roughly 50 — 150 kpc) vs average gas temperature for cool-core clusters (crosses) and non- 
cool-core clusters (diamonds), which are re-calculated based on the data of Zhang et al. (2007) and Sanderson et al. (2009). The coiTesponding best-fit power 
law models given by Sanderson et al. (2009) are also shown as dashed and dotted lines, respectively. The entropy-temperature distribution of the off-centre 
cool gas clumps (Table 3) is also plotted. 
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